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With a semiclassical quasi-static model we achieve an insight into the complex dynamics of two 
correlated electrons under the combined influence of a two-center Coulomb potential and an intense 
laser field. The model calculation is able to reproduce experimental data of nitrogen molecules for a 
wide range of laser intensities from tunnelling to over-the-barrier regime, and predicts a significant 
alignment effect on the ratio of double over single ion yield. The classical trajectory analysis allows 
to unveil sub-cycle molecular double ionization dynamics. 
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Within the strong-field physics community, there has 
been increasing interest on double ionization (DI) of 
molecules in intense laser pulses and a large variety of 
novel phenomena has emerged. The diatomic molecules 
show a much higher double ionization yield than the 
prediction of the single-active-electron (SAE) model by 
many orders of magnitude [l|, |2j, and DI yield as well 
as ionized-electron momentum distribution exhibit a 
strongdependence on molecular structure and alignment 
[3> EL TH- Experimental data indicate that a rescatter- 
ing mechanism is responsible for nonsequential double 
ionization (NSDI) caused by strong correlation between 
two electrons. In this process, which has been exten- 
sively investigated for atoms Q , an electron freed by tun- 
nelling ionization is driven by the laser electric field into 
a recollision with its parent ion. This essentially classi- 
cal rescattering picture implicitly suggests that the well 
synchronized recolliding electron burst with respect to 
the laser field be an alternative attosecond pulse source 
for the probe of molecular dynamics . However the de- 
tails of this pivotal recolliding event for molecules remain 
unknown. 

The complexity of the dynamics of the two correlated 
electrons responding to a two-center nuclear attraction 
and the laser force, on the other hand, poses a great 
challenge to any theoretical treatment. For instance, a 
time-dependent, three-dimensional quantum mechanical 
computation from first principles has not yet been ac- 
complished even for the simpler case of atoms Q. This 
leaves approximate approaches developed recently, such 
as one-dimensional quantum model [9(, many-body S- 
matrix and simplified classical methods How- 
ever, the complex electron dynamics which is crucial for 
molecular DI is still not fully explored and the theoretical 
results can not account for experimental data quantita- 
tively. In this letter, we employ a feasible semiclassical 
theory, i.e. an ab initio 3D calculation including classical 
rescattering and quantum tunnelling effects, providing an 
intuitive way of understanding the complex dynamics in- 
volved in the molecular DI. Our calculation is capable of 
reproducing unusual excess DI rate for a wide range of 



laser intensities quantitatively (see Fig. 1), thus consol- 
idating the classical rescattering view of molecular DI. 
In particular, with classical trajectory analysis, we are 
able to unveil the sub-cycle dynamics behind molecular 
DI and predict a significant influence of molecular align- 
ment on the ratio of double over single ion yield. 

The model we propose here is in the spirit of that of 
semiclassical treatment of DI of atoms in high-intensity 
field [12] • We consider a molecule composed of two nu- 
cleus and two valence electrons interacting with an in- 
frared laser pulse. When the laser intensity is smaller 
than a threshold value (see Fig. []}, one electron is re- 
leased at the outer edge of the suppressed Coulomb po- 
tential through quantum tunnelling (Fig. Ufa)) with a 
rate ru(to) given by molecular ADK formula fl3]. The 
initial position of the tunnelled electron can be de- 



10" 



10" 1 r 



10" r 



10" J r 



10' 



«* Exp. Data from C.Cornaggia 
— ADK theory 
-□ — Our model 




Threshold intensity 



0.1 1 
Peak laser intensity (1 15 W/cm 2 ) 



FIG. 1: (color online). Comparison between DI dataJH] and 
theory for nitrogen molecule. 0.185 PW/cm 2 is the threshold 
intensity separates the tunnelling and over-the-barrier regime 
as schematically plotted in Fig. 2. In the calculation, the laser 
frequency u) of 0. 05695a. u. and the number of optical cycle of 
37 are chosen to match the experiments of Cornaggia To 
our knowledge, the results from our calculations are the first 
to be in good agreement with experimental data for a wide 
range of laser intensities from tunnelling to over-the-barrier 
regime. 
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FIG. 2: (color online). Schematic representation of the ion- 
ization mechanisms addressed in this work, (a) Tunnelling 
ionization, (b) Over-the-barrier ionization. 



rived from the equation, — — — — + f \ — - — h-dr + 

n ' r a i r bl J n— r | 

J pl — zie{to) — with x\ = y\ = 0. The wave- 
function \l/ is given by the linear combination of the 
atomic orbital-molecular orbital (LCAO-MO) approx- 
imation [T3|. The initial velocity of tunnelled elec- 
tron is set to be (v± cos ip, v± sincp, 0), where vj_ is the 
quantum-mechanical transverse velocity distribution sat- 

isfying w(v±)dv± = P e( [ o) - exp( ' y £ ^ g > )dv±, 

and is the polar angle of the transverse velocity uni- 
formly distributed in the interval [0, 2tt] [12|. For the 
bound electron, the initial position and momentum are 
depicted by single-electron microcanonical distribution 
(SMD) H3,F(r 2 ,p2) = kS[I p2 - p 2 2 /2 - W(r a2 , r b2 )}, 
where k is the normalization factor, I P 2 denotes the 
ionization energy of molecular ions, and W(r a 2,ri,2) — 
— 1 /r a 2 — 1/ rb2 is the total interaction potential between 
the bound electron and two nuclei. 

The above scheme is only applicable when the laser 
intensity is lower than the threshold value [14] ]. To 
give a complete description of the DI of molecu- 
lar system for the whole range of the laser inten- 
sities (see Fig. 1), one need to extend the above 
model to the over-the-barrier regime (Fig. 2b). 
This is done by constructing the initial conditions 
with double-electron microcanonical distribution (DMD) 
QJ], i.e., F(ri,r 2 ,pi,p 2 ) = 2[/a(ri,Pi)/a(r2,p 2 ) + 

2 

//3(ri,p 1 )/a(r 2 ,P 2 )], with f aJj (r,p) = kS[I pl - \ - 
W(r a ,r b ) - V a , p (r)}, where ^(r) =^[1 - (1 + 

Kn, i a)e -2Kr6 ' a ] represents the mean interaction between 
the electrons, k can be obtained by a variational calcula- 
tion of the ionization energy of molecules. 

The subsequent evolution of the two-electron system 
with the above initial conditions is simulated by the 

,2 

classical Newtonian equations of motion: = e(t) — 
V(V^ e + V ee ). Here index i denotes the two different 
electrons, V^ e and V ee are Coulomb interaction between 
nuclei and electrons and between two electrons, respec- 
tively. V* e = - ^T,Vee = | ri l r2 | , where r al and 
rbi are distances between the ith electron and nucleus a 
and b. The above Newtonian equations of motion are 
solved using the 4-5th Runge-Kutta algorithm and DI 
events are identified by energy criterion. In our calcula- 



tions, more than 10 5 weighted (i.e., by rate ro(io)) clas- 
sical two-electron trajectories are traced and a few thou- 
sands or more of DI events are collected for statistics. 
Convergency of the numerical results is further tested by 
increasing the number of launched trajectories twice. 

The above model is applied to the study on DI of 
molecular nitrogen. We first calculate the ratio be- 
tween the double and single ionization yield with re- 
spect to the peak laser intensities from 5 x 10 13 W/cm 2 
to 1 x 10 15 W/cm 2 . In Fig. Q] the calculated results is 
compared with that of recent experiments of Cornaggia 
et.al. 0] and a good agreement is obtained for such a 
broad range of laser intensities. 

In the following, we proceed to explore the correlated 
electron dynamics responsible for DI of molecular nitro- 
gen. The classical trajectory method allows us to se- 
lect out the individual DI trajectories and back analyze 
their dynamics in detail. The typical electron trajecto- 
ries are shown in Fig. 3, presented in an energy ver- 
sus time plot. The threshold value of 0.185 PW/cm 2 
separates the DI data into two parts. When the peak 
laser intensity is below this value, there exist two dom- 
inant processes responsible for emitting both electrons, 
namely, collision-ionization(CI) and collision-excitation- 
ionization(CEI), as shown in Fig. EUa)(b), respectively. 
For CI, the tunnelled electron is driven back by the oscil- 
lating laser field to collide with the bound electron near 
its parent ion causing an instant (~ attosecond) ioniza- 
tion. For CEI, DI event is created by recollision with 
electron impact excitation followed by a time-delayed (~ 
a few optical periods) field ionization of the excited state. 
When the laser intensity is above the threshold value, 
over-the-barrier ionization emerges. In this regime we 
observe more complicated trajectories for DI processes. 
Except for CI (Fig. 3(c)) and CEI (Fig. 3(d)) trajecto- 
ries similar to tunnelling case, there are multiple-collision 
trajectories as shown in Fig. 3(e), (f) as well as collision- 
less trajectory of Fig. 3(g). In Fig. 3(e) and (f), initially 
two valence electrons entangle each other, experience a 
multiple-collision and then emit. The four types of tra- 
jectories shown in Fig. 3(c-f) represent the dominant pro- 
cesses of DI in the plateau regime from 0.185PW/cm 2 to 
0.5PW/cm 2 , each of them accompanied by one or multi- 
ple times of collisions between two electrons [l7|, llM [l9| . 
However, above 0.5PW/cm 2 , DI is dominated by a colli- 
sionless sequential ionization whose typical trajectory is 
represented by Fig. 3(g). In this regime results from our 
model agree with ADK theory. 

The analysis of trajectories of electron-electron pairs 
may provide insight into the complicated dynamics of DI 
with sub-cycle time resolution, and the important infor- 
mation is revealed by the laser field phase at the mo- 
ments of collision and ionization [2(| [H|. We choose 
three typical laser intensities, 0.12PW/cm 2 , 0.4PW/cm 2 
and lPW/cm 2 , representing the tunnelling, plateau and 
sequential ionization regime, respectively. 

Fig. 0Ja) shows the diagram of DI yield versus laser 
phase at the moment of closest collision. In the tun- 
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FIG. 3: (color online). Typical energy evolution of the elec- 
tron pair in different laser intensity regime. Vertical dashed 
lines indicate the moment when collision between electrons 
emerge. 




FIG. 4: (color online). DI yield vs laser phase when (a) 
the two electrons become closest; (b) both the electrons are 
ionized, at different laser intensity O.I2PW/cm 2 (triangle), 
0.4PW/cm 2 (circle) and lPW/cm 2 (square), respectively, (c) 
The relationship between the correlated momentum and the 
delay time at O.I2PW/cm 2 (upper), 0.4PW/cm 2 (mid) and 
f.OPW/cm 2 (low), respectively. 



nelling regime (i.e., 0.12PW/cm 2 ), we note that the col- 
lision can occur throughout most of the laser cycle and 
the peak emerges slightly before the zeroes of the laser 
field. This is consistent with the prediction of simple- 
man model I22II and recent results from purely classical 
calculation [23| . However, for the other two cases, the col- 
lision between the two correlated electrons occurs mainly 
at peak laser field. This is because the ionization mecha- 
nism changes at the transition to over-the-barrier regime, 
where both electrons rotate around the nuclei and their 
distance could be very close before one of them is driven 
away by the external field. 

Fig. H^b) plots DI yield as a function of the laser phase 
at the instant of ionization. Most DI occurs around the 
maximum of laser field for different intensity regime. In- 
terestingly, for the tunnelling case, we observe a peak 
shift of ~ 30° off the field maximum. With assuming 
that the colliding electron leaves the atom with no signif- 
icant energy and electron-electron momentum exchange 
in final state is negligible [24|, the parallel momentum 
kj 2 of each electron results exclusively from the accel- 
eration in the optical field: fc[ 2 = ±2\fU p sin uti on [2lj |. 
The above shifted peak indicates the accumulation of the 
emitted electrons at kj = fc 2 = ±VU p in the first and 

third quadrants of parallel momentum plane (fej , fcl ). It 
is consistent with the experimental data of Ref. [{J (see 
their Fig. 2). 

Fig. IHc) shows the phase angle of momentum vector 



(fci, fe 2 ) with respect to the delayed time between the 
closest collision and ionization. The integration over the 
phase angle gives total DI yield versus the delayed time. 
In all three cases we observe a long-tail up to several op- 
tical periods. For the sequential ionization of lPW/cm 2 , 
it means that the second electron is slowly (i.e., waiting 
for a few optical cycles) ionized after the first electron is 
deprived from nuclei by the laser field. In the tunnelling 
regime, the long-tail indicates that CEI mechanism is 
very pronounced for the molecular DI and contributes to 
- 80% of the total DI yield. 

This observation is different from purely classical sim- 
ulation [23[, where CI effect is believed to be overesti- 
mated. Our results, however, are consistent with experi- 
mental data for Ar atom [2(| , where ionization potential 
and laser field parameters are close to our case. The rea- 
son is stated as follows: For the intensity of 0.12PW/cm 2 , 
the maximal kinetic energy of the returned electron is 
3.17£/p = 0.85a.u., still smaller than the ionization en- 
ergy of Nq~. Even with the assistance of the Coulomb 
focusing [23, it is not easy for the returned electrons to 
induce too many CI events. Furthermore, such time delay 
might provide more physics beyond simple rescattering 
scenario. Recently a statistical thermalization model has 
been proposed for the nonsequential multiple ionization 
of atoms in the tunneling regime [26] . This model shows 
that sharing of excess energy between the tunnelled elec- 
tron and the bound electrons takes some time, resulting 
in a time delay on attosecond time scale between recolli- 
sion and ionization. Our simulation upholds this picture 



4 




30 60 

Molecular Alignment 



30 60 
Molecular Alignment 



90 



FIG. 5: (color online). The molecular alignment dependence 
of DI ratios for laser intensity of 0.15PW/cm 2 . 



of attosecond electron thermalization: on upper panel of 
Fig. Q|c) , two bright spots are observed at a similar time 
delay on subfemtosecond time scale for CI trajectory. 

The regular patterns in upper and lower panels of Fig. 
|4fc) exhibit that the ejection of electrons in the same- 
hemisphere and opposite-hemisphere emerge alternately 
with respect to the delayed time. For a time delay of 
odd half-lascr-cycles, two electrons emit in the same di- 
rection. In contrast, they emit in the opposite direction 
for an even half-laser-cycles time delay. In mid of Fig. 
Etc), on the other hand, the irregular pattern emerges 
as the signature of complicated multiple-collision trajec- 
tories for DI in the plateau regime. It implies that the 
trajectories of two electrons entangle with each other be- 
fore DI ionization occurs and the electrons' motion might 
be chaotic [27(. 

When the light intensity is high enough, it has been 
consensus that DI behavior of atoms is determined by 
essentially electron physics in the presence of laser field 
fl2L l23l | . Good correspondence between our theoretical 



calculations and experimental data confirms the valid- 
ity of the above picture in molecular DI case. In our 
model, after tunnelling electrons travel much of the time 
in the intense laser field like a classical object and solely 
electron collision physics determines the fate of DI of 
molecules. However, the inherent nuclear degree of free- 
dom of molecule do manifest themselves as the significant 
alignment effect in our model. To clearly demonstrate it, 
we calculate the ratios between double and single ioniza- 
tion at different molecular alignment angles. Main results 
are presented in Fig. [H It shows that, i) The ratio be- 
tween DI and single-ionization yield is less for perpendic- 
ular molecules than that of parallel molecules; ii) This 
anisotropy becomes more dramatic for a shorter laser 
pulse. Further explorations show that molecular align- 
ment also significantly affects the correlated momentum 
distribution of emitted electrons. Details will be pre- 
sented elsewhere p8l |. 

In summary, we exploit a semiclassical quasi-static 
model to achieve insight into the correlated electron dy- 
namics in molecular DI under the relevant experimental 
conditions, i.e., highly nonperturbative fields with fem- 
tosecond or shorter time resolution. Our calculation un- 
veils sub-cycle dynamics behind molecular DI and pre- 
dicts a significant influence of the molecular alignment on 
the ratio of double over single ion yield. Because molec- 
ular alignment is controllable with present technique Q 
the above results can be regarded as our theoretical pre- 
diction which may be tested in future experiments. 
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